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ABSTRACT 
A  major  hindrance  to  the  use  of  algebraic  techniques  for  the 
design  of  linear  control  systems  is  the  non-linearity  of  the  equations 
describing  the  system.  The  equations  for  a  general  class  of  control 
systems  are  shown  to  be  linear  in  certain  groupings  of  variables 
contained  in,  or  added  to,  the  systems.  A  design  technique  based 
on  the  solution  of  these  linear  equations  is  outlined,  and  a 
digital  computer  program  which  will  accomplish  the  design  is  pre- 
sented. 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA  93S43-5101  Table  of  contents 

Section  Page 

1.  Introduction  11 

2.  Algebraic  Methods  13 
Method  of  Mitrovic'  14 
Coefficient  Matching  15 

3.  Algebraic  Solution  of  a  General  Problem  18 
Steady-State  Error  Specification  24 
Bandwidth  Specification  25 
Multiple  Specifications  26 

4.  Limitations  of  the  Method  29 
Steady-State  Error  for  a  Ramp  Input  29 
Multi-path  Feedback  30 

5.  Sample  Solutions  32 

6.  Conclusions  75 
APPENDIX  Program  COMPTECH  77 


List  of  Tables 
Table  Page 

1.  Matrix  System  of  Equations  for 

Compensated  Feedback  Control  System.  23 


List  of  Illustrations 

Figure  Page 

1.  £  and  a)  in  the  s -Plane  13 

2.  Basic  Feedback  Control  System  18 

3.  Basic  Feedback  Control  System  with 

Cascade  Compensation  20 

4.  Control  System  with  Two  Feedback  Paths  30 

5.  Block  Diagram  for  Example  5.1  32 

6.  Block  Diagram  for  Example  5.2  41 

7.  Block  Diagram  for  Example  5.3  52 

8.  Block  Diagram  for  Example  5.4  63 


Table  of  Symbols  and  Abbreviations 

s  a  complex  variable  (»  =»  o  +  .w) 

X(s)  a  transfer  function;  a  ratio  of  polynomials  in  the 

variable  s. 

Xn(s)  the  numerator  of  X(s) 

Xd(s)  the  denominator  of  X(s) 

Xn,  the  coefficient  of  the  kth  power  of  s  in  Xn(s) 

f(x,y,z)  a  function  of  the  variables  x,  y  and  z 

Kt  variable  forward  gain 

T(s)  forward  transfer  function 

K  variable  feedback  gain 

H(s)  feedback  transfer  function 

C(s)  filter  compensation  transfer  function 

Ess  steady -state  error 

G„  gain  in  decibels  specified  in  a  bandwidth  requirement 

uD 

w  frequency  specified  in  a  bandwidth  requirement 

C  damping  ratio,  cosine  of  an  angle  in  the  s-plane 

w  undamped  natural  frequency,  a  radial  distance  in  the 
s-plane 


1.   Introduction 

Linear  control  system  design  may  be  generally  described  as  the 
synthesis  of  physical  components  or  systems,  active  or  passive,  whose 
dynamic  response  will  satisfy  specified  criteria.  Present  design 
procedure  is  essentially  mathematical,  since  the  response  of  dynamic 
systems  is  conveniently  described  by  a  set  of  differential  equations. 
The  problem  is  then  to: 

1)  formulate  differential  equations  whose  solution  will 
describe  the  desired  response, 

2)  relate  coefficients  of  the  equations  to  the  parameters  of 
a  physical  system, 

3)  obtain  a  solution  of  the  differential  equations  in  terms 
of  these  parameters. 

Production  of  systems  designed  in  this  manner  is  frequently  economi- 
cally unfeasible,  since  components  required  may  be  expensive  to  man- 
ufacture in  limited  numbers.  A  more  practical  approach  is  to  utilize 
readily  available  components,  combining  them  and  adjusting  their 
variable  parameters  to  obtain  a  system  with  the  desired  dynamic 
response. 

Analysis  of  systems  thus  synthesized  may  be  accomplished  in  the 
time  domain  by  direct  solution  of  the  system  differential  equations, 
but  the  more  usual  approach  is  to  perform  analysis  in  the  complex 
frequency  domain  (s-plane),  by  means  of  the  Laplace  Transformation. 
System  equations  thus  transformed  may  be  approached  with  several  well- 
known  analysis  techniques,  notably  the  Evan's  root  locus  procedure, 
or  in  the  frequency  domain  (s  ■  ju),  the  Bode,  Nyquist  or  Nichols 
methods.  Application  of  these  techniques  to  the  design  procedure  is 
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restricted,  however,  since  they  are  not  easily  described  analytically, 
they  are  essentially  graphical  and  therefore  time-consuming,  and  the 
approach  to  design  is  largely  trial  and  error.  All  suffer  the  further 
disadvantage  that  changes  in  system  response  may  be  predicted  on  the 
basis  of  variation  of  only  one  parameter. 

Several  recent  techniques,  to  be  discussed  in  further  detail, 
have  succeeded  in  overcoming  the  more  obvious  difficulties  listed 
above.  The  transformed  differential  equations  describing  the  system 
and  the  specifications  prescribing  system  response  are  reduced  to  a 
set  of  algebraic  functions  of  the  variable  parameters.  The  algebraic 
relationships  produced,  however,  are  non-linear,  and  for  complex 
systems  may  prove  difficult  or  impossible  to  solve.  This  paper  out- 
lines a  numerical  procedure  for  the  solution  of  these  non-linear 
equations,  and  presents  a  digital  computer  program  which  obtains 
solutions  for  a  general  class  of  control  device  problems. 
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2.  Algebraic  Methods 

Application  of  algebraic  techniques  requires  the  reduction  of 
the  transformed  system  equations  to  a  "transfer  function,"  usually 
the  ratio  of  two  polynomials  in  the  variable  s,  which  relates  the 
output  of  the  system  to  an  arbitrary  input.  The  roots  of  the  denom- 
inator polynomial,  known  as  the  "characteristic  equation,"  determine 
the  stability  and  response  time  of  the  system,  factors  of  principal 
concern  to  the  designer.  For  the  system  to  be  stable,  i.e.  finite 
response  to  a  finite  input,  it  is  necessary  that  all  roots  of  the 
characteristic  equation  have  negative  real  parts.  It  can  be  shown 
that  the  time  response  of  the  system  is  determined  by  the  values  of 
the  roots,  and  is  frequently  dominated  by  the  root  or  pair  of  roots 
whose  real  parts  have  the  least  negative  value,  i.e.  when  plotted  in 
the  s-plane,  the  roots  located  closest  to  the  origin.  The  dominant 
roots  are  designated  in  polar  coordinates  by  their  radial  distance 
from  the  origin  (cotj),  and  the  cosine  (5)  of  the  angle  formed  by  the 
second  quadrant  radius  vector  and  the  negative  real  axis. 

s-plane 

s  «■  a  +  .  to 

J 


E,     s  cos9 


Figure  1 


£  and  to  in  the  s -Plane 
n 
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Method  of  Mitrovic' 

The  characteristic  equation  is  a  polynomial  of  the  form 

F(s)  -  F  sn  +  F   .  &n~l  + +  F0s2  +  F-sV  F        (2.1) 

n     n-1  2      1     o      v   ' 

or  in  summation  notation 

n     k 
F(s)   =  E  F  s  (2.2) 

k=0  K 
where  the  polynomial  coefficients  (F,  's)  are  functions  of  the  system 
parameters.  Mitrovic' [1]  has  shown,  through  a  transform  of  variable 
and  s-plane  mapping  procedure,  that  two  coefficients  of  F(s)  may  be 
stated  as  functions  of  the  remaining  coefficients  and  of  the  %   and 
wn  of  a  pair  of  specified  roots.  The  result  of  this  procedure  is  two 
equations  of  the  form 

S   (-l)k  Fkwnk  d,   (O  -  0  (2.3) 

k=0       k     k_1 

I        (-l)k  Fk  oo/  0,k(O  -  0  (2.4) 

k=0 

where  <j>,  -,(?)  and  4>ip(^)  are  recurring  Chebishev  functions,  which  for 

a  specified  £  may  be  precomputed  from  the  relationship 

W(C>  "2?*k(?)    +    *k-l(e)    =    °  (2-5) 

with  <j>    (C)     -     0 
o 

and     <J>    (S)     =     -1 

If  numerical  values  are  assigned  to  all  but  two  of  the  coeffi- 
cients (F,  )  and  to  £  and  co   equations  (2.3)  and  (2.4)  are  reduced 
to  a  pair  of  simultaneous  equations  in  two  unknown  variables,  i.e. 
the  coefficients  which  have  not  been  specified.  Mitrovic'  originally 

chose  the  coefficients  F  and  F,  to  be  variables,  arriving  at  the 

o      1 

relations 
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Fo  =  fl  (Fk'  *»    ***>  <2'6> 

Fl  *  f2  (Fk'  ^    "ti  (2*7) 

Later,  Sikjak  generalized  the  Mitrovic1  method  to  show  that  a  solution 
can  be  obtained  for  any  pair  of  coefficients  [2] .  Thus,  a  relationship 
between  the  coefficients  of  the  characteristic  equation,  which  are 
functions  of  the  system  parameters,  and  the  location  of  a  desired  pair 
of  roots  is  established.  Designating  the  variable  system  parameters 
as  (a..,  a«,  — —  a  ),  we  have  for  any  two  coefficients 

Fx  =  f3(al»  a2»  """  az)     "     f4(Fk»  *>'  wn)  (2,8) 

Fy  =  f5(al'  a2»  —  az)  =  f6(Fk»  ?•  °°n>  <2-9> 

The  process  may  be  repeated  to  specify  additional  pairs  of  roots, 
thus  establishing  additional  relationships  of  the  form  of  (2.8)  and 
(2.9).   If,  however,  not  all  roots  of  the  characteristic  equation 
are  specified,  the  method  provides  no  information  on  the  location  of 
undesignated  roots,  and  does  not  guarantee  either  stability  or  the 
dominance  of  a  desired  pair  of  roots. 

Coefficient  Matching 

If  all  of  the  roots  of  F(s)  are  designated,  the  characteristic 
equation  may  be  stated  as 


F(s)  «   7r  ■<•  +  rk)  (2.10) 

k=»l 

where  (r. ,  r»,  r  )  are  the  roots  of  F(s).  The  resultant  coeffi- 

v  1*   2*      n' 

cients  of  the  expanded  polynomial  are  functions  of  the  roots. 

n     ,.    ,  v_n-l, 


F 


(s)  =  fn(r1,r2,-.-rn)s   +  fa-1  (r^r,.— «^i   +  — 

+  f2(r1.ra,-.rn)il2+£1<r1,rjl,T-Ti)«  +  f^r^r^-^) 

(2.11) 
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where  f,  (r..,r?,  - —  r  )  is  the  sum  of  the  products  of  all  possible 
combinations  of  (r.,^,  —  r  )  taken  k  at  a  time. 
For  example: 

Fk(s)  =  fk  (ri,r2,  — rn) 

-  (r1r2r3— rk)  +  (r^r^-r^)  +  jr^..^) 

+  Crir3r4™rk+1?  +  ---  (2'12) 

In  summation  notation 

x     n 

*V(S>  -2   [  t.  k  CO]  ;  x  =  n! (2.13) 


7AS)     -  Z   [  tt  k  (r.)]  ;  x  =  nl 

i=l   j=i    J         k!  (n-k).' 


Since  all  roots  are  designated,  the  F  ,'s  are  fixed  constants.  The 
result  is  two  expressions  for  F(s) 

n  i 

F(s)  ■  l      f  (a. ,a_— —  a  )s      Parameters  (2.14) 

i=0  1  L     l         2 

n 

F(s)  =  I      g  (r  ,r r  )s      Roots  (2.15) 

=0 

In  order  that  these  expressions  be  identically  equal,  coefficients 
of  the  same  powers  of  s  must  be  equal. 

Vk  =   fk(al>a2»—  az)sk  =  8k(rl'r2»— V^        (2'16) 
The  resultant  system  is  composed  of  n  simultaneous  equations  in  z  un- 
knowns. Complete  solution  of  the  system  of  equations  requires  that 
the  number  of  unknown  quantities  (z),  the  variable  parameters,  be  equal 
to  the  number  of  equations  (n),  the  order  of  the  characteristic  equa- 
tion.  If  insufficient  variables  are  available,  additional  parameters 
must  be  inserted  in  the  problem  in  order  that  this  requirement  be 
satisfied. 

The  principal  advantage  of  designation  of  all  roots  is  the  guar- 
antee of  stability  and  dominance  of  desired  roots.  By  designating  a 
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pair  of  roots  whose  values  give  the  desired  time  response,  and  loca- 
ting the  remaining  roots  in  such  a  way  as  to  insure  dominance  of  the 
original  pair  and  stability  of  the  system,  the  designer  completely 
specifies  the  time  response  of  the  system.  The  principal  disadvant- 
age is  the  necessity  for  addition  of  arbitrary  parameters  sometimes 
required  to  permit  solution  of  the  simultaneous  equations. 

A  disadvantage  common  to  all  algebraic  methods  is  the  inherent 
non-linearity  of  the  resultant  equations.  The  expressions  usually 
contain  products  and/or  powers  of  the  variable  parameters,  and 
assistance  of  a  digital  computer  is  required  for  solution  of  all 
but  the  most  simple  cases. 
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3.  Algebraic  Solution  of  a  General  Problem 

Given  a  dynamic  device  whose  transformed  time  domain  equation 
may  be  expressed  as  the  product  of  a  variable  gain  (Kt)  and  a  trans- 
fer function  T(s),  and  a  feedback  device  expressible  as  a  variable 
feedback  gain  (Kf)  and  the  transfer  function  H(s),  the  transfer 
function  F(s)  of  the  "closed  loop'1  control  system  comprising  the  two 
components,  as  arranged  in  Figure  2,  is 


Output  (s) 
Input  (s) 


F(s) 


Kt  T(s) 


1  +  Kt  T(s)  Kf  H(s) 


(3.1) 


Input(s) 


Figure  2 
Basic  Feedback  Control  System 


Both  T(s)  and  H(s)  are  the  ratio  of  polynomials  in  S. 

n.   ,   \  Tn  £a 

T(s) 


a-1  2 

Tn  8  +  Tn  .s   A  + +  Tn„s  +  Tn.,s  +  Tn 

-  Tp(s)  _  &       a~l  2       1      o 

=  Td(s) 


Td1.sb  +  Tdv  ,&h~l   + +  Td.s2  +  Td,s  +  Td 

D  D-l  Z  1        O 


H(s)  = 


Hd(s) 


(3.2) 


and  F(s) 


Kt  Tn(s) 
Td(s) 


1  +  Kt  Tn(s)  Kf  Rn(s) 
Td(s)  Hd(s) 


F(s)  reduces  to  the  ratio  of  two  polynomials 


(3.3) 
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F(s)  =  ISifii  =  Kt  Tn(s)-Hd(s) 

^  ;     Fd(s)     Td(s)-Hd(s)  +  Kt  Kf  Tn(s) -Hn(s)  K*'*} 

Forming  the  product  of  polynomials  where  indicated  defines: 

Tn  Hd(s)   =   Tn(s)«Hd(s) 

Td  Hd(s)   -   Td(s)«Hd(s) 

Tn  Hn(s)   -   Tn(s)»Hn(s) 
The  characteristic  equation  of  the  system  -  Fd(s)  is  the  sum  of  the 
polynomials  TdHd(s)  and  KtKfTnHn(s). 
Fd(s)   -  TdHd(s)  +  KtKfTnHn(s) 

-  TdHd  s*     +  TdHd   .s*"1  +  -—  +  TdHd.s  +  TdHd     (3.5) 
X  x-1  1         o 

+  KtKf(TnHn  «y  +  TnHn   ,sy"  + +  TnHn.,6  +  TnHn 

N    y        y-1  1        o 

In  most  devices  used  for  control  systems,  particularly  higher  order 

systems,  the  order  of  the  numerator  polynomial  is  much  less  than  the 

order  of  the  denominator.  Feedback  devices  seldom  have  numerator 

order  more  than  two  greater  than  denominator  order,  i.e.  second 

derivative  feedback.  Thus  the  order  of  the  characteristic  equation, 

and  the  number  of  roots  to  be  designated,  is  determined  by  the  order 

of  the  TdHd(s)  product.  By  expanding  the  designated  roots  into  a 

polynomial  R(s),  and  equating  coefficients,  the  system  of  simultaneous 

equations  for  Fd(s)  becomes 

(TdHdx  )  sX  =   Rx  sX 

(TdHd  .)  a*'1  =   R  .Bx"1 

v    x-1'  x-1 


(TdHd  +  KtKfTnHn   )»y  -   R  Sy 

y  y  y 

(TdHdy-1  +     KtKfTnHny^sy-1  -       K.^"1 

» 

(TdHd2  +     KtKfTnHn2     )s2  -       R^" 
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(TdHt^    -  KtKfTnHn1  )s 

(TdHd     =  KtKfTnHn   ) 
o  o 


R.  8 


R 


(3.6) 


If  all  polynomials  are  rationalized, 

TdHd       R    =   1 
x       x 

and  the  system  consists  of  x  equations  in  the  two  unknown  parameters 
Kf  and  Kt.  In  order  to  obtain  an  algebraic  solution,  an  additional 
(x-2)  variables  must  be  inserted  in  the  system.  The  addition  of 
"cascade  compensation"  in  the  form  of  a  multi-section  filter  net- 
work permits  introduction  of  the  necessary  mathematical  variables 
and  is  also  an  economical  method  of  inserting  variable  parameters 
in  the  physical  plant. 

Cascade  Compensation 


Input (s) 

— ^ 

v 

Kt 

Vfa. 

\ 

C(s) 

f 

-) 

JL\ei  / 

Output 

Kf 

IK«^ 

Figure  3 
Basic  Feedback  Control  System  with  Cascade  Compensation 
The  transformed  output/input  equations  C(s)  of  a  resistance- 
capacitance  filter  network  are  expressible  as  the  ratio  of  two 
equal  order  polynomials,  the  order  being  equal  to  the  number  of 


filter  sections  employed. 
C(s) 


Cn(s) 
Cd(s) 


Cn  ■"  +  Cn  ,sm"  + +  Cn„s  +  Cn,s  +  Cn 

m      m-1  2       1      o 


(3.7) 
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The  transfer  function  of  the  compensated  closed  loop  system,  for 
the  configuration  shown  in  Figure  3,  becomes 

Input  (s)        r^s;       1  +  Kt  T(s).C(s)  Kf  H(s)  KJ'0) 

which  reduces  to 

ws)  =  Ms)   =  Kt  TnHd(s)»Cn(s) 

v  '     Fd(s)     TdHd(s)*Cd(s)  +  KtKf  TnHn(s).Cn(s)     K***} 

The  order  of  the  characteristic  equation  is  now  determined  by  the 

TdHd(s)»Cd(s)  product,  thus  the  addition  of  each  filter  stage 

increases  by  one;  the  order  of  Fd(s),  the  number  of  simultaneous 

equations  to  be  solved,  and  the  number  of  independent  variables 

required.   If  the  order  of  TdHd(s)  is  (x),  and  (m)  filter  sections 

are  used,  the  order  of  Fd(s)  will  be  (x  +  m),  The  addition  of  each 

filter  section  (i)  introduces  two  independent  variables  into  the 

system,  Cn.  and  Cd . .  Thus  (m)  filter  sections  introduce  (2m) 

variables. 

For  solution  of  the  system  of  equations,  the  product  of  the 

independent  parameters  (KtKf)  may  be  considered  as  a  single 

algebraic  variable.  Addition  of  (m)  filter  sections  to  a  system 

whose  TdHd(s)  is  of  order  (x) ,  produces  (x  +  m)  equations  in 

(2m  +  1)  variables.  Equating  the  number  of  variables  and  number 

of  equations 

x  +  m  ■  2m  +  1 

m  -  x  -  1  (3.10) 

For  the  system  given  in  Figure  3,  the  number  of  filter  sections 

required  is  determined  by  the  product  of  the  denominator  polynomials 

of  the  forward  and  feedback  path  transfer  functions. 
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The  system  of  simultaneous  equations  for  the  compensated  system 
Is  obtained  by  designating  the  required  (x  +  m)  roots,  expanding  the 
roots  to  form  the  polynomial  R(s),  forming  the  TdHd(s)«Cd(s)  and 
KtKfTnHn(s)»Cn(s)  products,  and  equating  coefficients.   Table  1 
illustrates,  in  matrix  form,  the  equations  for  a  system  having  a 
third  order  Td(s)»Hd(s)  product,  a  first  order  Tn(s)»Hn(s)  product, 
and  the  required  second  order  filter.  The  format  of  the  matrix 
equations  is  the  same  for  higher  order  systems  if  the  requirements 
that 

1)  The  number  of  filter  sections  used  be  one  less  than  the 
order  of  Td(s)-Hd(s), 
and  2)  The  order  of  Tn(s)«Hn(s)  be  less  than  the  order  of 

Td(s)«Hd(s) 
are  fulfilled. 

Note  that  the  equations  for  systems  of  this  type  are  linear  in 
certain  individual  parameters  (Cd  ,  Cd.  —  Cd  ),  and  in  groupings 
of  the  other  parameters  (KtKfCn  ,  KtKfCn.,  —  KtKfCn^.  Thus  the 
system  of  equations  may  be  solved  by  linear  methods  for  these 
individual  and  grouped  parameters. 

Again,  if  all  polynomials  are  considered  rationalized, 

TdHd  Cd   «=  R  ,    =1 
x   m      x+m 

so  Cd  =1.   If  Cn  is  defined  as  equal  to  1,  the  group  variable 
m  m 

KtKfCn  becomes  simply  the  product  KtKf.  Having  solved  the  linear 
m 

system  and  obtained  the  coefficients  of  the  denominator  polynomial 

(Cd  ,  Cd.,  Cd  .)  of  C(s)  directly,  the  coefficients  of  the 

numerator  polynomial  may  be  obtained  by  dividing  the  grouped  variables 

(KtKfCn  ,  KtKfCn,,  KtKfCn  ,),  by  the  value  of  KtKfCn  ,  i.e. 

v      o*       1*  m-1"   '  m* 
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by  KtKf.  With  the  numerator  and  denominator  polynomials 
available,  the  filter  design  becomes  a  problem  in  passive  network 
synthesis. 

A  solution  has  thus  been  obtained  for  all  physical  parameters 
except  the  individual  gains 9  Kt  and  Kf.  An  additional  specification 
is  required  in  order  to  obtain  an  independent  solution  for  values 
of  the  gains.  The  additional  criteria  considered  are  steady  state 
error  and  bandwidth. 

Steady  -  State   Error  (Ess) 

The  error  of  a  system  is  defined  as  the  difference  between 
an  input  command  signal  and  the  system  response. 

Error  (s)  ■  Input  (s)  -  Output  (s)  (3.11) 

From  Equation  3.8 

Output  (s)  -  Kt  T(s)-C(s) 

1  +  Kt  T(s)«C(s)  Kf  H(s)   inpuC  lS; 

The  error  as  a  function  of  the  input  signal  becomes 

_     ,  N      1  +  Kt  T(s)»C(s)fKf  H(s)  °  11   _     /  \  /*  io\ 
Err°r  (S)  =  1  +  Kt  T(»)-C(a)  Kf  H(s)  '   Injmt  (s)  (3*12) 

which  reduces  to 

Error  (s)  _   TdHd(s)»Cd(s)  +   KtKf  TnHn(s)*Cn(s)  °  KtTnHd(s)«Cn(s) 
Input  (s)       TdHd(s).Cd(s)  +  KtKf  TnHn(s)«Cn(s)  ,3  13) 

For  a  unit  step  displacement  input 9  Input  (s)  is  1/s.  Error  (s) 

for  a  step  input  becomes 

/   \  1  i  Kt  TnHd(s)»Cn(;s)  fr.   ... 

Error  (s;     g    1  -  TdHd(s).Cd(s)  +  KtKf  TnHn(s).Cn(s)J  <"*«1£^ 

Steady-state  error  is  defined  as  the  error  of  the  system  after  all 
transient  effects  have  taken  place,  ©r  mathematically,  as  time 
expired  after  application  of  the  input  signal  approaches  infinity. 
Application  of  the  Final  Value  Theorem  to  Equation  3.14  gives 
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Ess  =  11m  s  Error  (s) 
s-*  o 


Ess  =  lim 
s-*o 


r     Kt  TnHd(s)»Cn(s) ~~| 

|_     TdHd(s)»Cd(s)  +  KtKf  TnHn(s) •Cn(s) J 


_        TdHd  •  Cd  +  TnHn  KtKfCn  -  TnHd  •  KtCn 

Ess  =       o   o       o o o o_  (~   -c\ 

TdHd#Cd   +  TnHn  •KtKfCn  l-»«"J 

o   o        o       o 

Note  that  steady-state  error  as  expressed  above  is  defined  for  unity 

input.  Error  is  relative  to  input,  and  maximum  steady-state  error 

desired  in  the  system  is  specified  as  a  fraction  of  the  input  signal. 

Steady-state  error  is  found  to  be  a  function  of  two  variables 

previously  defined,  Cd  and  KtKfCn  ,  and  an  additional  variable, 

KtCn  .  Values  for  Cd  ,  KtKfCn  and  KtKf  having  been  previously 

obtained,  the  solution  for  Kt  Cn  is 

■  o 

(1  -  Ess)  [TdHd  #Cd   +  TnHn  .KtKfCn  1 

KtCno     =     TnHd (3'16) 

o 

KtKfCn 

and  Kf  -  -KtcT^     *   Kt  -  W  (3-17) 

o 
Values  of  all  unknown  parameters  have  been  obtained,  and  a  system 
is  synthesized  which  will  meet  specified  requirements  for  stability, 
root  dominance,  and  steady-state  error. 

Bandwidth  (BW) 

The  bandwidth  of  a  system  is  defined  as  the  frequency  range 
over  which  the  response  of  a  control  system  is  within  acceptable 
limits  [1],   It  is  frequently  given  as  that  frequency  (w)  at  which 
the  magnitudes  of  the  output  response  and  the  input  signal  have  a 
specified  ratio  (G),  usually  expressed  in  decibels  (db). 

If  F(s)  is  the  closed  loop  transfer  function  of  a  control 
system,  bandwidth  may  be  found  from 
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/output   (luiM 

Gdb     "     20  lo*10  [  input     (jM)J~     20  lc^IK^)]                        (3.18) 

where   oo  is  the  desired  bandwidth. 

a  2 

Fn(  .  oj)            Fn  (  .  ui)     +  —  +  Fn   (  .  u)  +  Fn  (  .  u)  +  Fn 

G     ■     F(.o))     ■     ! — 5*-r     a     — *  <1 — -i=--J LJ 2 

j    '              Fd(ju)             Fdb(j'^)b  +  —  +  Fd2(jco)2+  Fdx(    a))   +  FdQ 


Fn(aQreal     +     .1  Fn(oi): 
Fd(w)real     +     j  Fd(w)imag. 


2 

oj)real)       + 

2  2 

oj)real)       +     (Fd(w)imag.) 


(3.19) 


The  bandwidth  expression  is  seen  to  be  highly  non-linear,  however, 
this  obstacle  is  not  insurmountable.  Recalling  that 

„,   N  Kt  TnHd(.u)*Cn(.(») 

F(,w)   -   vj      \i 

J        TdHd (  .  u) •  Cd (  .  u)  +  TnHn(  .  w)#  KtKfCn(  .  u) 
J      J  J         J 

and  that  values  for  the  variables  Cn,8  Gd„,  and  KtKfCn.  are  available 

i    i  i 

from  the  solution  of  the  simultaneous  equations,  Kt  may  be  obtained 
from  the  bandwidth  criterion.  If  gain  (G)  is  specified  at  frequency 


Kt 


CTdHd  (  .  w)»  Cd (    a)     +     TnHn(    <$•  KtKfCn(  .  <^ 
TnHd(.(^.Cn(„w)  J  <^-^ 


and         Kf     ■     KtKf  (3.21) 

Kt 

Again,  values  of  all  unknown  parameters  are  obtained,  and  the  system 

design  will  meet  stability,  root  dominance  and  bandwidth  specifications, 

Multiple  Specifications 

Either  a  steady-state  error  or  a  bandwidth  specification  is 

sufficient  in  itself  to  permit  solution  for  all  variable  parameters 

and  design  of  the  system.  However,  if  both  error  and  bandwidth 
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are  specified,  solution  for  one  does  not  guarantee  satisfaction  of 
the  other  requirement.  A  slightly  different  approach  must  be  taken. 

No  requirement  has  been  placed  on  the  location  of  the  non- 
dominant,  or  non-specified  roots,  other  than  that  they  have  negative 
real  parts,  and  that  they  be  sufficiently  removed  from  the  origin  of 
the  s-plane  to  insure  dominance  of  the  desired  roots.  This  gives  the 
designer  a  great  deal  of  latitude  in  the  choice  of  extraneous  roots 
which  are  added  to  the  problem  as  a  result  of  cascade  compensation. 

The  system  of  linear  equations  used  in  the  solution  for  para- 
meters is  obtained  by  equating  coefficients  of  the  characteristic 
equation.  On  one  side  of  each  equation  is  an  algebraic  function 
of  the  variable  parameters;  on  the  other  side  is  the  numerical 
coefficient  found  by  expanding  the  roots  to  form  the  polynomial 
R(s).  Changing  the  location  of  one  or  more  roots  changes  the  values 
of  the  variable  parameters.  Parameter  values  obtained  are  those 
which  will  place  the  roots  of  the  characteristic  equation  exactly 
where  specified.  The  steady-state  error  specification  is  satisfied 
by  determining  the  relationship  of  the  forward  gain  (Kfc)  to  other 
system  parameters.  The  bandwidth  is  found  from  the  equation 


G 


Kt  TnHd(  w)  Cn(.u))  ,^  ^\ 


RCjto) 


Thus  the  bandwidth  of  the  system  is  a  function  of  the  location  of 
the  roots,  but  the  steady-state  error  solution  is  not.  By  judicious 
choice  of  his  root  locations,  the  designer  may  be  able  to  meet  both 
error  and  bandwidth  specifications,  or  to  determine  whether  a  band- 
width specification  conflicts  with  a  stability  or  root  dominance 
criterion. 
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The  process  of  finding  the  proper  root  locations  is  iterative, 
laborious,  and  time-consuming 9   requiring  a  complete  solution  of  the 
system  equations  after  each  movement  of  the  roots.  For  higher  order 
systems,  the  assistance  of  a  digital  computer  may  be  necessary. 
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4.  Limitations  of  the  Method 

Steady-State  Error  for  a  Ramp  Input 
Steady-state  error,  as  defined  in  Section  3,  is  the  limit  as 
s  approaches  zero  of  the  product  of  s  and  the  error  function 
Error  (s),  where 

Errorfs)  =  Iuimt(3)[*dHd(')'Cd(8)  +  KtKfTnHn(s)*Cn(s)  -  KtTnHd(s)>Cn(s)1 
arrows;   inpuc^s;^        TdHd(s).Cd(s)  +  KtKf  TnHn(s)Cn(s)        J 

For  a  unit  ramp  input,  i.e.  constant  velocity  signal  of  unity 

2 
magnitude,  Input  (s)  ■  1/s  .  Steady-state  error  becomes 

s  fTdHd(s)»Cd(s)  +  KtKfTnHn(s)*Cn(s)  -  KtTnHd(s)»Cn(s)  ~| 
Ess  =siim  ji   |^         TdHd(s).Cd(s)  +  KtKf  TnHn(s)Cn(s)      J 

This  limit  approaches  infinity  as  s  approaches  zero,  unless  an  s  can 

be  factored  from  the  numerator  polynomial  of  the  term  in  brackets,  i.e. 

.  .  TdHd  ♦  Cd   +  TnHn  •  KtKfCn   -  TnHd  -  KtCn   =0        (4.2) 
oo       oo        oo 

Equation  4.2  may  be  set  equal  to  zero  in  several  ways: 

1)  The  possibility  of  requiring  that  Cn  or  Cd  be  zero  must 
be  disregarded,  as  this  would  fix  the  value  of  a  variable 
required  for  solution  of  the  system  of  simultaneous  equa- 
tions. 

2)  Hn  may  be  required  to  be  zero,  i.e.  tachometer  feedback, 
but  the  steady-state  error  expression  resulting  from  this 
procedure  introduces  two  new  variables  KtCn  and  KtCn. 

in  a  single  equation 

Ess  -  f(Cd  ,  Cd,,  KtKxCn  ,  KtCn  ,  KtCn.) 
v  o*   1'       o'     o*     1' 

and  solution  for  independent  Kt  is  impossible. 
Kt  Cn  Kt  Cn. 

Kt  "  -sr2     *    Kt  ■  -zr1 

o  1 
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3)  Td  may  be  required  to  be  zero,  i.e.  a  pole  at  the  origin 

of  the  open  loop  transfer  function.  Then 

Tn  Hn  KtKfCn   -  Tn  Hd  KtCn   ■  0 
o  o      o      o  o     o 


and 


Kf 


Hd    :  i.e.  unity  feedback. 
o 

Hn 


If  this  procedure  is  used  K£  is  no  longer  an  independent 

variable  and  solution  is  impossible. 

4)  If  Hd  is  required  to  be  zero, 

Tn  Hn  KtKfCn   =  0 
o  o      o 

or  Tn  must  also  be  zero,  and  an  s  factors  from  tta®  daaomin* 

ator  as  well  as  from  the  numerator  of  Error  (s). 

To  date,  the  author  has  found  no  expression  for  steady-state 

error  to  a  ramp  input  which  is  compatible  with  the  system  of 

simultaneous  equations  previously  derived. 

Multi-path  Feedback 


Input (s) 


Output (s) 


Figure  4 
Control  System  with  Two  Feedback  Paths 
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The  transfer  function  for  the  control  system  illustrated  in 
Figure  4  may  be  written 

*KS)  1  +  Kt  T(s)-C(s)  [Kf  H(s)  +  Kg]  KH'J) 

which  reduces  to 

F(S)  =  Kt  Tn  Hd(s)»Cn(s) ~  .. 

K   '  TdHd(s)*Cd(s)  +  KtKf  TnHn(s).Cn(s)  +  KtKg  TnHd(s)»Cn(s)  v  *  J 

The  system  of  simultaneous  equations  obtained  by  equating  coefficients 

of  the  characteristic  equation  will  be  linear  in  the  variables  (Cd  , 

Cd  ,  -—  Cd  ),  (KtKfCno,  KtKfCn  ,  —  KtKf)  and  (KtKgCnQ,  KtKgCn^  — 

KtKg) .  The  equations  may  be  solved  to  obtain  values  for  the 

variables  composed  of  grouped  parameters.  However,  for  this  system, 

and  other  multi-path  feedback  systems  of  various  configurations,  it 

has  not  proved  possible  to  solve  for  independent  values  of  the 

independent  parameters*   For  example,  in  the  system  above  KtKfCn  , 

KtKf,  KtKgCn  ,  and  KtKg  are  independent  variables.  But  the  value  of 

Cn  obtained  from 
o 

Cn   ■  KtKfCn 
0  o 

KtKf 

is  not  necessarily  the  value  of  Cn  obtained  from 

o 

Cn   -  KtKgCn 
o       ^  o 

KtKg 
Systems  of  this  type  have  defied  solution  by  this  method. 
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5.  Sample  Solutions 
Example  5.1 


*V©-  Kt 


(s+l)  (s-2)  (t+8)  (s+15) 


0(8) 


Figure  5.1 
Block  Diagram  for  Example  5.1 

The  open  loop  transfer  function  given  above  is  unstable  due 
to  the  pole  in  the  right  half  plane  (s  ■  2).  Design  a  filter 
compensator  and  find  values  for  forward  and  feedback  gain  which 
will  make  the  system  shown  in  Figure  5.1  absolutely  stable. 

Dominant  roots  are  specified  by  £  =  0.707  and  go  ■  2.5. 

No  additional  steady-state  error  or  bandwidth  requirements 
are  specified. 
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INPU1  DATA 
#»*#«♦»*#♦## 

OPEN  LOOP  TRANSFER  FUNCTION  -  T(S) 

NUMERATOR  GF  OPEN  LOOP 
TRANSFER  FUNCTION  -  TN(S) 

TN(S)=   l.OOOOE+OQ 

DENOMINATOR  OF  OPEN  LOOP 
TRANSFER  FUNCTION  -  TD(S) 

NUMBER  OF  ROCTS  =  ^ 

ROOT    MO.                     12  3  4 

REAL    PT.  -.llL+ul          .20E+01  -.80E+G1  -.15E+02 

IMAG    PT.          ..COE  +  CO          .JOF  +  ::1  .CO. .  +  00  .COE+00 

COEFFICIENTS    Of     EXPANDED    POLYNOMIAl 

td(S)=  i.:~cof+:o  S(  4) 

+  2.20C0E+  1  S(  3) 

+  9.50CGE+01  S(  2) 

+  -1.6600E  +  92  S(  1  ) 

+- 2.400GF+  2 

NOTE  -  PROGRAM  HAS  READ  POLE(S) 
IN  THE  RIGHT  HALF  OF  THE 
COMPLEX  PLANE. 
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FEECBACK  PATH  TRANSFER  FUNCTION  -  H(S) 
(0)  ORDER  FEEDBACK  IS  SPECIFIED. 

NUMERATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HN(S) 

HN(S)=   l.OOOOE+00 

DENOMINATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HD(S) 

HD(S)=   l.OOOOE+CO 


CLOSED  LOOP  TRANSFER  FUNCTION 
G(S)  -  UNCOMPENSATED 

NUMERATOR  CF  CLOSED  LOOP 
TRANSFER  FUNCTION  -  GN(S) 

GN(S)=  KT*   1.U000E+00 


DENOMINATOR  OF  CLOSED  LOOP 
TRANSFER  FUNCTION  -  GDIS) 

GD(S)=   1-.90GOE+C0  S(  4) 

+  2.200OE+01  S(  3) 

+  9.50C0L  +  01  S(  2) 
♦-1.6600E+02  S(  1) 
♦  (-2.AOOOE+02  ♦  KT«KX»  l.OOGOE+CG) 
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STEADY  STATE  FRRUR  SPEC  I  F  I  C  AT  I  ON 

NO  STEADY  STATr  ERROR  REQUIREMENT 
IS  SPECIFIED. 


BANDWIDTH  SPECIFICATION 

NO  BANDWIDTH  REQUIREMENT  IS  SPECIFIED. 

NOTE  -  DATA  INDICATES  NO  BANDWIDTH  OR  PARAMETER 

DEPENDENT  STFADY  STATE  ERROR  SPECIFICATIONS. 
UNDER  THESE  CONDITIONS  THE  PROGRAM  WILL  NOT 
SOLVE  FOR  FORWARD  GAIN  AND  FFEDBACK  GAIN 
SEPARATELY.  SOLUTION  IS  FOR  THE  PRODUCT  OF 
THE  GAINS.   (KT*KX  ) 


ROOTS  SPECIFIED  FOR 
CLOSEC  LOOP  SOLUTION 


ETA  =   .7C7CE+C 


NUMBER  (IF  ROOTS  =  2 


OMEGA  =   2.500  E+00 


POOT  NO. 


REAL  PT.   -.1CF  +  -  1   -.18E+  '•! 


IMAG  PT .    .18E+01   -.18E+~1 
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SOLUTION   OF    THE    PROBLEM 


PROGRAM    REQUIRES       7    ROOTS    FOR    SOLUTION. 


SOLUTION    OF    THE    SYSTEM    OF    LINEAR    EQUATIONS 


VARIABLE 


VARIABLE    (     1) 


VARIABLE    (    2) 


VARIABLE    (     3) 


VARIABLE    (     4) 


VARIABLE    (     5) 


VARIABLE    (    6) 


VARIABLE    (    7) 


VARIABLE    (     8) 


PARAMETERS 
CD(O) 

com 

CD(2) 

CD(3) 

KT*KX»CN(0) 

KT«KX»CN(1) 

KT«KX#CN<2) 

KT»KX»CN(3) 


RELATIVE  ERROR 

7.358E-09  PERCENT 

.OO0E-O1  PERCENT 

.OOOE-Ol  PERCENT 

-4.180E-09  PERCENT 

-1.032E-07  PERCENT 

3.466E-07  PERCENT 

1.393E-06  PERCENT 

2.421E-06  PERCENT 


NOTE  -  RELATIVE  ERROR  STATED  ABOVE 
MEASURES  THE  ACCURACY  OF  THE 
SOLUTION  FOR  EACH  VARIABLE. 
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SOLUTION  FOR  INDIVIDUAL  PARAMETERS 

NUMERATOR  QF  CASCADE  COMPENSATION 
CN(S)=   1.0OOOE+C0  SI  3) 

+  L.1867E+01  S(  2) 

♦  3.A32CE+0!  S(  1) 

+  2.7105E+01 

ZEROS  OF  THE  FILTER 

ROOT  NO.  1          2  3 

REAL  PT.  -.8QE+01  -.13E+01  -.25E+01 

I  MAG  PT.  .COil+flC  .OOE+QO  .OOE+OO 

NOTE  -  ACCURACY  OF  ROUT  FINDING  SUBROUTINE 
HAS  BEEN  CHECKED  BY  EXPANDING  THESE 
ROOTS,  AND  COMPARING  THE  COEFFICIENTS 
OBTAINEC  WITH  THE  COEFFICIENTS  OF  THE 
ORIGINAL  POLYNOMIAL. 

ERROR  NCTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROUTS  IS  CONSIDERED  ACCURATE. 
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DENOMINATOR   OF   CASCADE    COMPENSATION 
CDIS>*      l.OOOOE+OO    S(    31 

♦  1.5117E+01  S(  2) 

♦  1.4418F+02  S(  I) 

♦  2.9161E+02 

POLES  OF  THE  FILTER 

ROOT  NO.        1         2         3 

REAL  PT.   -.63E+01   -.63E+01  -.26E+01 
IMAG  |t;   -.85E+01    .85E*01    .OOE+OO 

NOTE  -  ACCURACY  HAS  BEEN  CHECKEOtAS  BEFORE. 
ERROR  NOTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR  ' 
ROOTS  IS  CONSIDERED  ACCURATE. 

PRODUCT  OF  GAINS  (KT»KX)  «  5.4508E+03 

FORWARD  GAIN  (KT)  IS   5.4508E*03f 
WITH  UNITY  FEEDBACK  GAIN.   (KX=1.0) 
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CLOSED  LOOP  TRANSFER  FUNCTION 
OF  SOLVED  PROBLEM 

NUMERATOR  POLYNOMIAL 

GN( S)=   5.4508E+03  S(  3) 

+  6.4686E+C4  S(  2) 

♦  1.8707E+C5  S(  1) 

+  1.A775E+C5 

ZEROS  OF  THE  TRANSFER  FUNCTION 

ROOT  NO.  1          2          3 

REAL  PT.  -.80E+01   -.13E+01   -.2SF+C1 

IMAG  PT.  .CCE+CO    .OCE+OO    .OOE+GO 

NOTE  -  ACCURACY  FAS  BEEN  CHECKED,AS  BEFORE. 
ERROR  NOTED  IS  LESS  THAN  CNE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 
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DENOMINATOR  POLYNOMIAL 

GD(S)=   l.OOCOE+00  S(  7) 

♦  3.7118H+01  S(  6) 
+  5.7177E+C2  S(  b) 
+  4.7338E+C3  S(  A) 

♦  2.2814E+u4  S(  3) 
+  6.4826E+Q4  S(  2  ) 
+  1.0406E+05  S(  1  ) 
+  7.7761E+C4 

POLES  OF  THE  SYSTEM  AND  ROOTS 
CF  THE  CHARACTERISTIC  EOUATICN 

ROOT  NO.        1  2  3  4  5 

REAL  PT.   -.18E+01   -.18E+01   -.53E+G1   -.53E+01   -.71E>01 

IMAG  PT.    . 18E+01   -.18E  +  '>1    .IBE+OO   -.18C+00    .18E+00 

ROOT  NO.        6  7 

REAL  PT.   -.71E+01   -.88E+01 
IMAG  PT.   -.18f+:0    .OOE+OO 
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Example  5.2 


Ks) 


+ 


=\ 

Kt 

s  +  5 

»       — 

J" 

s4  +  33s3  +  336s2  +  1880s  +  7200 

0(s) 

Kf 

s  (s  +  7) 

s2  +  30s  +  216 

Figure  5.2 
Block  Diagram  for  Example  5.2 

Cascade  compensation  and  the  tachometer  feedback  device  given 
are  to  be  used  in  the  closed  loop  system  shown  in  Figure  5.2.  The 
maximum  permissible  steady-state  error  for  a  step  input  is  to  be 
10  percent. 

Dominant  roots  are  specified  by  ^  »  0.5  and  oj  =6. 


El 


No  bandwidth  requirement  is  specified. 
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INPUT  DATA 

OPEN  LOOP  TRANSFER  FUNCTION  -  T(S) 

NUMERATOR  OF  OPEN  LOOP 
TRANSFER  FUNCTION  -  TN(S) 

TNCS)*   1.0000E400  S(  1) 

♦  5.OOCOE+O0 

DENOMINATOR  OF  OPEN  LOOP 
TRANSFER  FUNCTION  -  TO(S) 

TDIS)«  KOOOOE  +  OO  SI  4) 

♦  3.3000E+CI  S(  3) 
+  3.3600E+02  S(  2  1 

♦  1.8800E  +  03  S(  1  ) 

♦  7.2000H-03 

ROOTS   OF    TOTS) 

ROOT  NO,        1          2          3  4 

REAt  PT.   -.90E+01  -.20E+02  -.20E+01  -.20E+01 

IMAG  PT.    .OOE+OO  .OOE+OO  -.60E+01  .60E+01 
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FEECRACK  PATH  TRANSFER  FUNCTION  -  H(S) 
(1)  ORDER  FEEDBACK  IS  SPECIFIED. 

NUMERATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HN(S) 

FN(S)=   i.OOOOE+OO  S(  2) 

+  7.GOCQE+CC  S(  1  ) 

+   .OOOOE-Ol 

DENOMINATOR  OF  FFFDHACK  PATH 
TRANSFER  FUNCTION  -  HD(S) 

HD(S)=   l.GOCOe+OO  S(  2) 

+  3.00COF  +  C  1  S(  1  ) 

+  2.16COt+02 
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CLOSET,    LOUP    TRANSFER    FUNCTION 
G(S)     -    UNCOMPENSATED 

NUMERATOR    OF    CLOSED    LOOP 
TRANSFER    FUNCTION    -    GN(S) 

GN(S)=    KT*(        l.OOCPC+OO  S(  3) 

+    3.5  JO    t  +  Ol  S(  ? ) 

+    3.660OE+O2  S(  1) 

+    1.0800E+03  ) 


CENOMINATOR    OF    CLOSED    LOUP 
TRANSFER    FUNCTION    -    GDIS) 

GD( S)=  1. . jOOE+OG    S(    6) 

+  6.3GOCE+01    S(     5 ) 

+  1.5420F+Q3    S(    4) 

+  (     1.9088E  +  H4    +    KT»KX*    l.GOGOE+OOl     S(     3) 

+  (     1.3618E+05    +    KT«KX»    L200OE+rl)     S(    2) 

+  (     6.2208E  +  :5    +    KT*KX*    3.50  -OE+.-l)     S(     1) 

+  (     1.5552E06    +    KT*KX«       ,OOOOE-:il 
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STf-ADY  STATE  ERROR  SPEC  I  F  IC  AT  1  flN 


ERROR  SPECIFIED  (E^ROR  /  INPUT)  IS   .1COOE+UO 


BANDWIDTH  SPECIFICATION 


NO  BANDWIDTH  REQUIREMENT  IS  SPECIFIED, 


RO(  TS  SPECIFIED  FOR 
CLOSED  LOOP  SOLUTION 


ZETA    =       .5TQC L+r    • 


OMEGA    =       6.00C  E+OI) 


NUMBER  OF  ROOTS  =  2 


ROOT    NO. 


REAL    PT.       -.30F>"1       -.30E+G1 


IMAG    PT. 
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SOLUTION  OF  THE  PROBLEM 


PROGRAM  REQUIRES  11  ROOTS  FOR  SOLUTION- 


STEADY  STATE  ERROR  SOLUTION 

STEADY  STATE  ERROR  IS   •10O0E+OO 

RELATIVE  SPECIFICATION 
ERROR  IS   3.727E-06  PERCENT 

ERROR  SPECIFICATION  IS  SATISFIED. 
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SOLUTION  OF  THE  SYSTEM  OF  LINEAR  COUATIONS 


VARIABLE 


VARIABLE  (  1) 


VARIABLE  (  2) 


VARIABLE  (  3) 


VARIABLE  (  4) 


VARIABLE  (  b) 


VARIABLE!  (  6) 


VARIABLE  (  7) 


VARIABLE  {  8) 


VARIABLE  (  9) 


VARIABLE  (1.1) 


VARIABLE  (11) 


VARIABLE  (12) 


VARIABLE  (13) 


PARAMETERS 
CD(C  ) 
CO(  1) 
CD(2I 
CD(  3) 
CO(  A) 
CD(  5) 

KT*KX«CNO) 
KT*KX*CN(1) 
KT*KX»CM(2) 
KT*KX*CN( 3) 
KT»KX*CN(4) 
KT»KX»CN(5) 
KT»CN(  0) 


RELATIVE  ERROR 

.0C.'E-01  PERCENT 

-1.A61E-09  PERCENT 

-^.589E-09  PERCENT 

] .759E-09  PERCENT 

.OD'E-01  PERCENT 

.00?E-01  PERCENT 

-•5.42JE-09  PERCENT 

-1.27  E-Q8  PERCENT 

-3.431E-08  PERCENT 

9.962E-08  PERCENT 

P.723E-07  PERCENT 

3.H3E-06  PERCENT 

^.727E-06  PERCENT 


NOTE  -  RELATIVE  FRROR  STATED  ABOVE 
MEASURES  THE  ACCURACY  OF  THE 
SOLUTION  FOR  EACH  VARIABLE. 
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SOLUTION  FOR  INDIVIDUAL  PARAMETERS 

NUMERATOR  CF  CASCADE  COMPENSATION 

CN(S)=   l.nOCOE+OC  S(  5) 

♦  6.7668F+01  S(  4) 

♦  1.7775E+03  S(  3) 

♦  2.2639E+04  SI  2) 

♦  1.3997E+05  S<  1 > 
+  3.3682E+C5 

ZEROS  OF  THE  FILTER 

ROOT  NO.        1          2          3          A  5 

REAL  PT.   -.20L+?2   -.12E+02  -.18E+32  -.84E+01  -.92E+G1 

IMAG  PT.    .COE+JO    .OOE+OO  .OOE+GO    .OOE+OO  .OOE+GO 

NOTE  -  ACCURACY  CF  ROOT  FINDING  SUBROUTINE 
HAS  BEEN  CHECKED  BY  EXPANDING  THESE 
ROOTS,  AND  COMPARING  THE  COEFFICIENTS 
OBTAINED  WITH  THE  COEFFICIENTS  OF  THE 
ORIGINAL  POLYNOMIAL. 

ERROR  NOTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 
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DENOMINATOR    OF    CASCADE    COMPENSATION 

CD(S)=  l.OOGOE+OO  S(  5) 

+  7.2'Oje  +  Ol  S(  4) 

+  2.     591r>03  S(  3) 

+  2.7312E+'">4  S(  2) 

+  1.7124E+05  S(  1) 

+  4. 16341  + j5 

POLFS    OF  THE    FILTER 

ROOT    NO.                     1                           2  3                          4  5 

REAL    PT.  -.24L+.<2       -.24E+02  -.89E+C1  -.78E+01  -.78E+C1 

IMAG    PT.  -.HE  +  02         .11E+02  .OOE+uO  -.26E+91  .26L+J1 

NOTE    -    ACCURACY    FAS    BEEN    CHECKED, AS    BEFORE. 
ERROR    NOTED    IS    LESS    THAN    CNE    PERCENT 
FOR    ALL    COEFFICIENTS.     SOLUTION    FOR 
ROOTS     IS    CONSIDERED    ACCURATE. 

FORWARD    GAIN     (KT)     -     1.6020E+03 
FEEDBACK    GAIN    (KX)    =    1.2011F+C0 
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CLOSED  LOUP  TRANSFER  FUNCTION 
(  F  SOLVED  PROBLEM 

NUMERATOR  POLYNOMIAL 

GN(S)=   1.6020E+03  S(  3) 

+  1.6447F+C5  S(  7) 

♦  7.2280L+06  S(  6) 
+  1.7734F+C8  S(  5) 

♦  2.6528E+G9  S(  4) 
+  2.4737E+10  S(  3) 
+  1.4G12E+11  S(  2) 
+  4.3966E+11  S(  1 ) 
+  5.8275h+ll 

ZEROS  OF  THE  TRANSFER  FUNCTION 

ROOT  NO.  12          3          4  5 

REAL  PT.  -.20E+02   -.18E+*<2   -.18C+02  -.12E+02  -.12E+02 

IMAG  PT.  .OOE+OO   -.46E-J2    .46E-02  -.36E-02  .36E-02 

ROOT  NO.        ft  7  8 

REAL  PT.   -.84E+01   -.92E-K>1   -.50E+01 
IMAG  PT.    .OOE+OO    .OOE+O^    .OOE+GO 

NOTE  -  ACCURACY  HAS  BEEN  CHECKED, AS  BEFORE. 
ERROR  NOTFD  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 


50 


DENOMINATOR    I'ULYNOMAL 

GDC S)«  l.OCC :t  +  "0    SC 11  ) 

+  1.35C-OE+02    S(  10) 

+  6.1371E+03    SC     9) 

+  2.8907E+C5    SC    8) 

+  6.73*  ri[  +  T6    SC     7) 

♦  1  .  :.  H  1  C  r  +  ~  8  S  (  6  ) 
+  1.2275E+C9  SC  5 ) 
+  9.9345E+C9    S(    4) 

♦  5.^6631+10  SC  3) 
+  2.2290F+11  SC  2) 
+  5.4  7991+11  S(  1) 
+  6.4750H+11 

POLES    OF    THE    SYSTEM    AND    ROOTS 
CF    THE    CHARACTERISTIC    EQUATION 


ROOT    NO. 


1 


REAL  PT.   -.30E  +  01   -.3uE+;l   -.90E+01   -.9CE+01   -.12E+02 
IMAG  PT.    .52E+01   -.52E+31    .52E+CO   -.52E+00    .52E+00 


ROOT  NO. 


10 


REAL  PT.   -.12E  +  '..~   -.15E+j2   -.15E+02   -. 18E+02   -.18E+C2 


IMAG  PT.   -.t-2L+. 


.52E+00   -.52E+°0    .52E+30   -.52E+00 


ROOT  NO. 


11 


REAL  PT.   -.21E+02 


IMAG  PT.    .0OE+C3 
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Example  5.3 


Msl 


T&- 


+   5 


Figure  5.3 
Block  Diagram  for  Example  5.3 


The  gain  of  the  open  loop  device  shown  in  Figure  5.3  is  less 
than  -20  db  for  all  frequencies.  Design  a  filter  compensator  and 
find  values  of  forward  and  feedback  gain  which  will  provide  positive 
gain  at  low  frequencies  and  a  bandwidth  of  zero  db  at  30  radians  per 
second. 

Dominant  roots  are  specified  by  E,   s  0.5  and  oo  =  10. 

No  steady-state  error  requirement  is  specified. 
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rNPUT     DATA 


OPfcN    LCfl      TRANSFER    FUNCTION    -    T(S) 

NUMERATOR    OF     OPEN    LOOP 
TRANSFFR    FUN(  HON    -    TN(  S) 

TN(S)=       6.3    :j.  +  ,1 

DENOMINATOR    (  r    UPE  J    LOUP 
TRANSFLR     FUi  f  T ION    -    TD(  S) 


TD(S)=       1.1 


+   >0    o(     A  ) 


+     3.  *   •'         +     1    S(      U 


+    2.73CO    *     2    S{     2 ) 


+    7 .  A    f  u  L  +  ,  2    S  (     1  ) 


4-       r        7 .' >  '  +     I  ~> 


ROUTS    OF    T  T  (  S ) 


ROUT    ,\|U. 


REAL    PT.       -,12t+02       -.15F+Q2       -.14F+J1 


IMAG    PT. 


.  )OE+i 


.COE+OO 


-.25E+01 

.OOE+OO 


53 


FEEDBACK  PATH  TRANSFER  FUNCTION  -  H(S) 
(0)  ORDER  FEEDBACK  IS  SPECIFIED, 

NUMERATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HN(S) 

FN(S)=   l.OOOOE+00 

DENOMINATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HD(S) 

HD(S)=   l.OOCOE  +  00  S(  1  ) 

+  5.OOOOE+C0 


54 


CLOSED  LOOP  TRANSFER  FUNCTION 
G(S)  -  UNCOMPENSATED 

NUMERATOR  GF  CLOSED  LOOP 
TRANSFER  FUNCTION  -  GN(S) 

GN( S )=  KT*(   6.3C00E+01  S(  1 ) 

+  3.1503E+02  ) 

DENOMINATOR  VF    CLOSED  LOOP 
TRANSFER  FUNCTION  -  f,D(S) 

GD(  S)=  l.cr.jr+co    S(    5) 

+  3.5-COh+Cl    S(    4) 

+  4.23C0F+  '?    S(     3 ) 

+  2.1  .DO;  +    3    S(     ?) 

+  4.?760F  +  r>3    S(     1  ) 

+  (     2.e8jOE  +  33    +    KT*KX*    6. 3QOOE+' 1 ) 


55 


STEADY  STATE  fcRRC,<  SPECIFICAT 

NO  STEADY  STATE  ERROR  REQUIREMENT 
IS  SPECIFIED. 


BANDWIDTH  SPECIFICATION 

GAIN  OF    .GOOE-Ol  DB   AT  FREQUENCY 
OF   3.000E+O1  RAD/SEC  IS  SPECIFIED. 


ROOTS  SPECIFIED  FOR 
CLOSED  LOt.P  SOLUTION 

ZETA  =   .500CE+00        OMEGA  =   l.OOO.E+Ol 

NUMBER  OF  ROCTS  =  2 

ROOT  NO.        1  2 

REAL  PT.   -.50E+01   -.50E+01 

IMAG  PT.    .87E+01   -.87E+01 
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SULU1 ION  OF  THE  PROBLEM 


PROGRAM  RECUIRi  S   9  ROOTS  FOR  SOLUTION. 


BANDWIDTH  SOLUTION 


GAIN  IS    .OCffc-Ol  DB   AT  FREQUENCY  3."OOE+01  RAD/SEC 


RELATIVE  ERROR  IS    .OOCfc-Jl  PERCENT 


BANDWIDTH  SPEC  If-  ICU  ION  IS  SATISFIED 
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THE  SYSTEM  OF  LINtAR  EQUATIONS 


PARAMETERS 


RELATIVE  ERROR 


VAR  [ ACL   (  1  ) 


C  D  (  0  ) 


-2.355E-09  PERCENT 


VARIABLE  (  2) 


CD(  1) 


-5.282F-09  PERCENT 


VARIABLE  (  3) 


CD!  2) 


-7.561E-08  PERCENT 


VARIABLE  (  A) 


CLH  3) 


-3.551E-07  PERCENT 


VARIABLE  (  f5) 


CD(4) 


1.539E-06  PERCENT 


VARIABLE  i     b) 


KT«KX»CN( 0) 


7.654E-05  PERCENT 


VARIABLE  (  7) 


KT»KX*CN (1) 


l.'.68E-03  PERCENT 


VARIABLE  (  8) 


KT*KX*CN(2) 


8.962E-03  PERCENT 


VARIABLE  (  9) 


KT»KX*CN(3) 


r.28jE-02  PERCENT 


VARIABLE  (1j) 


KT*KX*CN(4) 


2.441E-01  PERCENT 


NOTE  -  RELATIVE  ERROR  STATED  ABOVE 
MEASURES  THE  ACCURACY  OF  THE 
SOLUTION  FOR  EACH  VARIABLE. 
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SHLOTIUN  FOR  INDIVIDUAL  PARAMETERS 


NUMERATOR  UF  CASCADE  COMPENSATION 


cn( s )=  i.  r  t  +  :  s<  4) 


+  3.ai>C7L-  +  'M  S(  3) 


+  5 .  'i  2  1  vr  +  •: -2    S  (  2  ) 


+  3 . 4  03  4  E  +  0  3  S  (  1  ) 


4  8. '34  95'  +0  3 


ZEROS  OF  THE  Fit  IE 


ROOT  Mi'!. 


PFAL  PT.   -.12E+N2 


IMAG  PT. 


-.1dE+  2 
.0 3E+?0 


-.6LE+01   -.61E+  1 


-.36E+01    .  36E-K.1 


NtJTt  -  ACCURACY  ff   ii'l  FINDING  SUBROUTINE 
HAS  BEEN  CHECKED  BY  EXPANDING  THESE 
ROOTS,  AND  COMPARING  THr  COEFFICIENTS 
OBTAINED  hI7H  THE  COEFFICIENTS  OF  THL 
ORIGINAL  POLYNOMIAL. 

ERRHR  NOTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROUTS  lb  CONSIDERED  ACCURATE. 
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DENOMINATOR  OF  CASCADE  COMPENSATION 
CD(S)=   1.0024L+00  St  4) 
+  1.25C0E+02  S(  3) 

♦  6.3543F+03  S(  2) 

♦  1.6789L+05  S(  1) 
+  2.4427E+06 

POLES  OF  THE  FILTER 

ROOT  NO.  1  2          3  A 

REAL  PT.  -.47E+02  -.47E+J2  -.15h02  -.15E+02 

IMAG  PT.  -.17F+02  .17E+02  -.27^  +  02  .27E+32 

NOTE  -  ACCURACY  HAS  BEEN  CHECKED, AS  BEFORE. 
ERROR  NCTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 

FORWARD  GAIN  (KT)  =  A.2411E+Q4 
FEEDBACK  GAIN  (KX)  =  7.1267E+00 
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CLOSED    LOUP    TRANSFER    FUNCTION 
I F    SOLVED    PROBLEM 

NUMERATOR    POLYNOMIAL 

CN(S)=       2.6719'.  +06    S(  5  ) 

♦     1.1625^+08    S (  A  ) 

+    1.961    r  + ..  n   s  (  3  j 


+     1.6336E+10    S(     2  ) 


+  6.s  jii:-+i:  s<   l ) 


+     1.1422E+11 


ZEROS    (F     FK-    TRANSFER    FUNCTION 


ROOT    NO. 


REAL    PT.       -.M'+.l       -.61E+01       -.15E+02 


[MAG    PT.       -.36E+01  .36E  +  ',1  .OCE  +  00 


-.12E+02 
,00€+00 


-.50t+01 
•09E+00 


NOTE    -    ACCURACV    HAS     BEEN    CHECKEO»AS    BEFORE. 

ERROR    NCTED     IS    GREATER    THAN    ONE    PERCENT 
FOR    AT    LEAST    ONE    COEFFICIENT.     SOLUTION 
FUR    ROOTS     IS    NOT    CONSIDERED    ACCURATE. 
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DENOMINATOR  POLYNOMIAL 

GD(S)=   l.COC0fc+O0  S(  9) 

♦  1.60COE+02  S(  8) 

♦  1.1152E+Q4  S(  7) 
+  4.4527E+C5  S(  6) 

♦  1.1274E+07  S(  5) 

♦  1.8947E+08  S(  4) 
+  2.1474E+C9  S(  3) 
+  1.62C1F+1Q  S(  2) 
+  7.5736E+10  S(  1) 
+  1.6983E+11 

POLES  OF  THE  SYSTEM  ANO  ROOTS 
CF  THE  CHARACTERISTIC  EQUATION 

ROOT  NO.        1  2  3  4  5 

REAL  PT.   -.50E+01   -.50E+01   -.15E+02   -.15E+02   -.20E+02 

IMAG  PT.    .87E+C1   -.87E+01    .87E+00   -.87E+00    .87E+C0 

ROOT  NO.        6  7  8  9 

REAL  PT.   -.20E+02   -.25E+02   -.25E+02   -.30E+O2 
IMAG  PT.   -.87E+0C    .87E+00   -.87E+00    .OOE+OO 
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INPUT  DATA 


OPEN  LOCP  TRANSFER  FUNCTION  -  T(S) 

NUMERATOR  OF  OPEN  LOOP 
TRANSFFR  FUNCTION  -  TN(S) 

NUMBER  OF  ROOTS  =  2 

ROOT  NO.        1  2 

REAL  PT.   -.60E+01   -.13E  +  '2 


IMAG  PT. 


•  00 E  +  Ot  .  0 C  E  +  j  .; 


COEFFICIENTS  OF  EXPANDED  POLYNOMIAL 


TN( S)=   1. 3CCOC+00  S(  2 ) 


+  1.9jC jt+01  S(  1  ) 


+  7. 80. TOE +01 
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DENOMINATOR  OF  UPEN  LOOP 
TRANSFCR  FUNCTION  -  TD(S) 

NUMBER  OF  ROOTS  =  <r. 

ROOT  NO.         1 


REAL  PT. 


-.5C-E+01   -.5CE+01   -.9CE+G1   -.15E+G2 


IMAG  PT. 


;+O0    .30E+01   -.30E+01 


-OE+OO    .80E+01 


ROOT  NO. 


RFAL  PT  .   -.  IV  +■  2 


IMAG  PI  .   -.80F+.>1 


COEFFICIENTS  Of   rX'ANOED  POLYNOMIAL 


td(  s  )=     l. :     :-f+co  s(   6  ) 


+  a  .  g :  c  o  ;•-  +  o  l   s  (    5 ) 


+    9.8  300E-.+    <"!    Sf     4  ) 


+    9.5170E+U3    St     3) 


+    4.5;i6L+C4    S(     2) 


+     8.&434:  +  '"4    S(     1  ) 


+        .?JCOL-Jl 
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FEEDBACK  PATH  TRANSFER  FUNCTION  -  H(S) 
(0)  ORDER  FEEDBACK  IS  SPECIFIED. 

NUMERATOR  OF  FtEOBACK  PATH 
TRANSFER  FUNCTION  -  HN(S) 

HN(S)=   l.OOOOE+00 

DENOMINATOR  OF  FEEDBACK  PATH 
TRANSFER  FUNCTION  -  HD(S) 

HD(S)=   l.OOOOE+00 
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CLOSFC    LOl;P     TRANSFER    FUNCTION 
G(S)    -    UNCOMPENSATED 

NUMERATOR    OF    CLOSED    LOOP 
TRANSFER    FUNCTION    -    GN(S) 

GN(S)=    KT«(       l.OOOJE+GO    S(    2) 

+    i.9C'Ov'E+01     S(    1  ) 

+    7.80C0E+Q1     ) 


DENOMINATOR    CF    CLOSED    LOOP 
TRANSFER    FUNCTION    -    GD(S) 

cd( s)=     i. x : oe+oo  s(   6 ) 

+  4.9     : 'uE  +  '.U    S(     5  ) 

+  9. 83!  CI  +.2    S(     4  ) 

+  S. 5 1701:  +03    S(     3  ) 

+  (     4.5016E+J4    +    KT*KX*    l.CC    DE+fO)     S(    2) 

+  (     8.P434E+,4    +    KT»KX»    1.9000E+01)     S'(     1) 

+  (       .CGOGE-91    +    KT*KX«    7.80CGE+01) 


67 


STEADY  STATE  ERROR  SPECIFICATION 
ERROR  SPECIFIED  (ERROR  /  INPUT)  IS   .10G0E+00 


BANDWIDTH  SPECIFICATION 

GAIN  OF  -3.0C0E+00  DB   AT  FREQUENCY 
OF   2.CO0E+O1  RAD/SEC  IS  SPECIFIED. 


ROOTS  SPECIFIED  FOR 
CLOSED  LOUP  SOLUTION 

ZETA  =   .707CE+00        OMEGA  =   4.0000E+00 

NUMBER  OF  ROOTS  =  2 

ROOT  NO.        1  2 

REAL  PT.   -.28E+01   -.28E+C1 

IMAG  PT.    .28L+01   -.28E+01 


tt 


SOLUTION  UF  THE  PROBLEM 


*•»**«***»•*•«*•••■»«*•*•• 


PROGRAM  REQUIRES  11  ROOTS  FOR  SOLUTION 


BANDWIDTH  SOLUTION 


ITERATIVE  PROCEDURE  FOR  BANDWIDTH  SOLUTION 


PROCESS  IS  CONVERGING 


ITERATION 


GAIN-DB 


RELATIVE  ERROR 


31 


-2.64E+C0 


1.20E+C1  PERCFNT 
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-3.-2L+00    -6.65E-J1  PERCENT 


CAIN  IS  -3.C'2Qfc  +  O0  DB   AT  FREQUENCY  2.00E+G1  RAD/SEC 


RELATIVE  ERROR  IS  -6.655E-01  PERCENT 


BANDWIDTH  SPECIFICATION  IS  SATISFIED 


STEADY  STATE  ERROR  SOLUTION 

STEADY  STATE  ERROR  IS   .  lOO.E  +  '-M 

RELATIVE  SPECIFICATION 
ERROR  IS   3.274E-08  PERCENT 

ERROR  SPECIFICATION  IS  SATISFIED 
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SOLUTION  OF  THE  SYSTEM  OF  LINEAR  EQUATIONS 


VARIABLE 


VARIABLE  (  1) 


VARIABLE  (  2) 


VARIABLE  (  3) 


VARIABLE  (  4) 


VARIABLE  {  5) 


VARIABLE  (  6) 


PARAMETERS 

CD(O) 

CD(  1) 

com 

CD(3) 

CD(4) 

CD(  5) 

RELATIVE  ERROR 


-2.290E-09  PERCENT 


.OOOE-Ol  PERCENT 


.uOOE-Ol  PERCENT 


1.689E-09  PERCENT 


.OOOE-Ol  PERCENT 


3.595E-09  PERCENT 


VARIABLE  (  7) 


VARIABLE  (  8) 


VARIABLE  (  9) 


VARIABLE  (10) 


VARIABLE  (11) 


VARIABLE  (12) 


VARIABLE  (13) 


KT*KX*CN(0) 
KT*KX*CN(1) 
KT»KX»CN(2) 
KT*KX»CN<3) 
KT*KX»CN(4) 
KT*KX»CN(5) 
KT»CN(  0) 


1.032E-08  PERCENT 


1.033E-07  PERCENT 


1.555E-06  PERCENT 


1.347E-05  PERCENT 


8.095E-05  PERCENT 


A.292E-04  PERCENT 


-5.123E+09  PERCENT 


NOTE  -  RELATIVE  ERROR  STATED  ABOVE 
MEASURES  THE  ACCURACY  OF  THE 
SOLUTION  FOR  EACH  VARIABLE. 
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SOLUTION  FOR  INDIVIDUAL  PARAMETERS 

NUMERATOR  OF  CASCADE  COMPENSATION 

CN( S)=   1. 30COE+O0  S(  5) 

'♦  4.5730E+G1  S(  A) 

+  8.3424C+02  S(  3) 

+  6  .  9  9  7  9  E  +  0  3  S  (  2  ) 

+  2.7415t+)4  S(  1 ) 

+  4.7798E+  4 

ZcPOS  OF  THE  FfLTEr, 

ROOT  NO.         1           2  3          A          5 

REAL  PT.   -.15F  +  C2   -.15E+w2  -.90EO1  -.34E+,'l  -.34E+01 

IMAG  PT.   -.76f+  1    .76E+01  .OOE+CO  -.27E+C1  .27L+01 

NOTE  -  ACCURACY  CF  ROOT  FINDING  SUBROUTINE 
HAS  BEEN  CHECKED  BY  EXPANDING  THESE 
ROOTS,  AND  COMPARING  THE  COEFFICIENTS 
OBTAINED  WITH  THE  COEFFICIENTS  OF  THE 
ORIGINAL  POLYNOMIAL. 

ERROR  NOTED  IS  LESS  THAN  CNF  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 
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DENOMINATOR  CF  CASCADE  COMPENSATION 

CD(S)=   l.OOOOE+OO  S(  5) 

+  1.7057E+O2  S(  41 

+  1.1647t+04  S(  3) 

+  4.C007E+05  S(  2) 

+  4.0045E+06  S(  1) 

+  5.8576E+06 

POLES  OF  THF  FILTER 

ROOT  NO.        1          2          3  4  5 

REAL  PT.   -.80E+02   -.38E+02  -.38E+02  -.18E+01  -.13E+02 

IMAG  PT.    .CGE+OO   -.43E+02  .43E+02  .OOE+OO  .OOE+OO 

NOTE  -  ACCURACY  HAS  BEEN  CHECKED t AS  BEFORE. 
ERROR  NOTED  IS  LESS  THAN  ONE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 

FORWARD  GAIN  (KT)  =  2.6980E+C6 
FEEDBACK  GAIN  <KX)  =  1.1111E+00 
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CLOSED    LOOP     TRANSFER    FUNCTION 
H     SOLVED    PROBLEM 

NUMERATOR    POLYNOMIAL 

GN(  S  )=       2.6980t  +  ')6    S(  7) 

+     1.7464!  +08    S(  6  ) 

♦    4.8.-55C  +  09    S(  5 ) 


-»     7. 12691+10    S(     4) 


+    6,:-853F+ll    S(     3) 


+    3  .  }  1 2  1 1  +  1  2    S(    2) 


+    8.241 7t+ 12    S(    1 ) 


+    1.0059E+13 


ZEROS    OF    THF    TRANSFER    FUNCTION 


ROOT    NO. 


REAL    PT.       -.13E+0?      -.15E+02       -.15F+02       -.60E+01      -.90E+01 

IMAG  PT.    .COL+&0   -.76E+01    .76E+01    .00E+3C    .OOE+00 

ROOT  NO.        6  7 

REAL  PT.   -.34E+G1   -.34E+J1 
IMAG  PT.   -.27-  +  '.  1    .27E+H 

NOTE  -  ACCURACY  HAS  BEEN  CHECKED, AS  BEFORE. 
ERROR  NOTED  IS  LESS  THAN  CNE  PERCENT 
FOR  ALL  COEFFICIENTS.  SOLUTION  FOR 
ROOTS  IS  CONSIDERED  ACCURATE. 
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DENOMINATOR  POLYNOMIAL 
GD<S>  =   LOQCOE+OO  S(ll) 
+  2.1957E+02  S<10> 

♦  2.0988E+G4  S<  9) 
+  1.1480E+06  S(  8) 
+  3.9723C+07  S(  7) 

♦  9.08COE+08  S(  6) 

♦  1.3910E+10  S(  5) 
+  1.4210E+11  S(  4) 
+  9.4753E+11  S(  3) 
+  3.9646E+12  S(  2) 
+  9.6743E+12  S(  1  ) 
+  1.1177E+13 

POLES  OE  THE  SYSTEM  AND  ROOTS 
CF  THE  CHARACTERISTIC  EQUATION 

ROOT  NO.        1  2  3  4  5 

REAL  PT.   -.28E+01   -.28E+Q1   -.85E+01   -.85E+01   -.17C+C2 

IMAG  PT.    .28E+01   -.28E+G1    .28E+00   -.28E+OG    .28L+00 

ROOT  NO.        6  7  8  9         10 

REAL  PT.   -.17E  +  02   -.26E+02   -.26E+02   -.34E+02   -.34E+02 
IMAG  PT.   -.28E+00    .28E+00   -.28E+03    .28E+00   -.28E+0Q 

ROOT  NO.       11 

REAL  PT.   -.43E+02 

IMAG  PT.    .OOE+00 
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6.   Conclusion 

The  method  presented  has  been  shown  to  be  fast  and  accurate, 
and  capable  of  designing  systems  which  will  meet  specifications 
for  root  dominance,  steady-state  error,  and  bandwidth. 

The  principal  disadvantage  of  the  method  presented,  exclusive 
of  the  restrictions  mentioned  in  Section  4,  is  the  large  number  of 
filter  sections  required,  particularly  for  higher  order  open  loop 
systems.  The  filter  sections  tend  to  further  complicate  the 
originally  complex  systems. 

The  possibility  of  reducing  the  number  of  filter  sections 
employed  in  the  final  design  is  worthy  of  further  study.  Approxi- 
mation of  the  characteristics  of  the  cascade  compensation  may  be 
accomplished  with  a  filter  of  lower  order,  but  analysis  of  the 
effects  of  such  approximation  on  system  specifications  is  beyond 
the  scope  of  this  paper. 
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